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FOREWORD 


Installation  of  an  additional  sensor  on  the  nose  of  a  missile  may  require  a  complete  re¬ 
design,  whereas  mounting  on  the  side,  away  from  the  nose,  may  be  relatively  easy.  However, 
use  of  a  side-mounted  sensor  would  require  that  the  missile’s  heading  error  be  sufficiently 
large,  greater  than  a  specified  angle,  in  order  that  the  sensor  ‘see’  its  target.  This  note  exam¬ 
ines  the  kinematics  of  flight,  constrained  in  this  way,  for  use  against  a  stationary  radiating 
target  on  the  ground.  If  the  missile  is  required  to  descend  onto  the  target  from  above,  for 
the  purpose  of  continuous  imaging  the  plan,  view  of  the  target,  an  additional  constraint  is 
needed.  Calculation  of  the  miss  distance  achievable  under  different  flight  conditions  is  done 
to  access  the  feasibihty  of  a  side-mounted  sensor  for  this  purpose. 
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and  Aegis  Programs.  This  note  has  been  reviewed  by  Dr.  C.  F.  Fennemore,  Technical  Lead, 
Target  Tracking  and  Signal  Processing;  Technology  Program  Manager,  Robin  Staton;  and 
F.  P.  Rucky,  Digital  Systems  Branch. 
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ABSTRACT 


Terminal  guidance  under  the  constraints  that  (i)  the  heading  error  must  be  greater  than  a 
specified  angle  (in  order  that  the  sensor  see  the  target),  and  (ii)  the  magnitude  of  the  lateral 
acceleration  is  limited,  is  studied.  Only  the  kinematic  features  of  flight  are  considered.  It 
is  demonstrated  that  a  planar  trajectory  is  a,s  good  as  any  for  the  purpose  of  hitting  a 
stationary  target  on  the  ground.  In  order  to  hit,  or  come  close  to  the  target,  the  final 
segment  of  terminal  flight  must  violate  the  heading-error  constraint,  during  which  time  the 
target  is  invisible. 

If  there  is  an  additional  condition  {Hi)  that  the  missile  must  remain  above  the  target 
{e.g.,  for  imaging)  during  terminal  descent,  the  optimum  trajectory  is  a  conic  spiral.  As  in 
the  previous  case,  there  is  a  final  segment  of  blind  flight.  The  miss  distance  depends  on  the 
constraints  and  other  parameters  of  the  problem,  as  well  as  on  the  sensor  measurement  errors. 
Monte  Carlo  simulation  was  used  to  estimate  the  miss  distance  under  various  conditions. 
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CHAPTER  1 
INTRODUCTION 

In  passive  homing,  it  is  sometimes  advantageous  to  place  the  sensor  on  the  side  of  the 
missile  or  warhead  (e.^.,  on  a  ballistic  warhead  where  placement  on  the  nose  would  cause  a 
difficult  heating  problem).  However,  a  passive  sensor  on  the  side  of  the  missile  may  prevent 
looking  in  the  direction  of  motion.  In  such  cases,  an  alternative  homing  technique  where  the 
missile  velocity  is  constrained  away  from  the  look  direction  is  required.  A  guidance  law  for 
the  problem  of  engaging  a  quasi-stationary  target  on  the  groxmd  by  a  vertically-descending 
side-looking  missile  is  needed. 

Aerodynamics  is  to  be  completely  neglected.  Thus,  there  is  no  missile  body,  and  therefore 
no  look  direction  relative  to  any  body  axis.  The  velocity  vector  essentially  takes  the  place  of 
the  body  axis  and  the  heading  error  takes  the  place  of  the  look  angle.  If  the  target  is  located 
at  the  origin,  and  r  and  f  denote  the  instantaneous  missile  position  and  velocity,  while  ri^ 
and  £0  represent  the  initial  values  of  these  vectors,  then  the  heading  error  <{>  is  determined 

by 

cos^  =  — =— =  (1.1) 

rv 

(see  Figure  1-1)  where  r  stands  for  |r|  and  v  stands  for  |r|.  The  placement  of  the  sensor 
will  make  it  impossible  to  ‘see’  the  target  if  the  look  angle  is  too  small.  This  constraint  is 
expressed  by  requiring  that  the  heading  error  <f>  never  be  less  than  a  given  critical  angle  (j)c- 

In  this  simplified  analysis  it  will  be  assumed  that  v  remains  constant.  This  is  somewhat 
realistic  if  the  missile  is  dropping  toward  the  target  in  an  almost  vertical  trajectory,  in  which 
case  V  will  be  approximately  determined  by  the  missile  terminal  speed.  The  terminal  speed 
decreases  with  increasing  air  density,  but  this  small  effect  is  neglected.  It  will  also  be  assumed 
that  the  missile  is  capable  of  only  a  limited  lateral  acceleration  to  be  denoted  by  Oc,  which 
may  depend  on  the  speed  v. 

It  is  likely  that  no  trajectory  will  be  able  to  actually  hit  the  target  under  the  reasonable 
conditions  outhned  below.  If  this  is  the  case,  the  problem  becomes  one  of  approaching 
“near  enough”  to  the  target  to  destroy  it  with  a  proximity  fuse,  or  to  hit  part  of  a  large 
extended  target.  It  may  be  necessary  to  specify  another  parameter.  Dm,  the  maximum 
allowed  miss  distance  whereby  trajectories  approaching  within  Dm  of  the  target  are  deemed 
to  be  successful. 
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Figure  3-1.  The  Heading  Error, 
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CHAPTER  2 

PLANAR  TRAJECTORIES 


This  chapter  considers  the  terminal  flight  phase  during  which  the  missile  approaches  the 
target  along  a  planar  trajectory.  The  flight  previous  to  the  terminal  phase,  in  which  the 
missile  comes  from  a  more  distant  launching  point  and  arrives  to  the  vicinity  of  the  target, 
is  not  of  interest  and  is  not  discussed.  This  preliminary  flight  phase  can  be  governed  by  any 
guidance  scheme  including  Proportional  Navigation  (PN),  and  is  here  referred  to  as  the  PN 
phase.  The  Heading  Error  (HE),  denoted  by  <f>,  must  always  exceed  a  prescribed  minimum 
value  <f)c  in  order  that  the  sensor  see  the  target,  as  explained  in  Chapter  1. 

In  the  approach  (or  PN)  phase,  the  HE  will  generally  decrease  as  the  velocity  vector  tends 
to  point  closer  and  closer  toward  the  target.  As  soon  as  the  Heading  Error  (HE)  becomes  as 
small  as  (l>c,  the  trajectory  will  transition  to  Constant  Heading  Error  (CHE).  The  reason  for 
this  is  that  increasing  the  HE,  although  allowed  by  the  constraint,  wiU  degrade  the  trajectory 
and  delay  the  arrival  time  in  the  deterministic  case  or  increase  the  miss  distance  in  stochastic 
cases. 

During  a  CHE  trajectory  the  lateral  acceleration  magnitude  (LA)  increases  without 
limit  as  the  turn  radius  decreases.  This  continual  decrease  in  the  turn  radius  is  illustrated 
in  Figures  2-l(a)  and  2-l(b).  Since  the  CHE  trajectory  wiU  not  approach  near  enough  to 
the  target  without  unlimited  LA,  there  is  a  transition  to  a  CLA  (constant  LA)  trajectory, 
which  is  a  circular  arc.  During  this  final  CLA  trajectory  the  HE  is  generally  too  small  for 
the  target  to  be  visible  to  the  sensor,  meaning  that  the  missile  is  flying  blind  during  the 
terminal  CLA  phase. 

Let  acLA  tlie  LA  required  for  the  missile  in  its  current  state  to  transition  to  a  CLA 
trajectory  and  hit  the  target.  The  logic  determining  the  type  of  trajectory  required  can  be 
expressed  by  the  following  algorithm: 

while  LA  <  ac/2 
while  (f>  >  (j)c 

use  PN  (or  any  appropriate  approach  guidance) 
when  <f)  =  (f>c 
use  CHE 
\i  <i>  <^c 

This  condition  should  ordinarily  never  be  encountered. 

end 

when  LA  =  ac/2 
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transition  to  CLA,  whereupon  LA  becomes  equal  to  Cc 

end 

TRAJECTORY  SHAPE 

It  is  shown  below  that  the  planar  trajectory  is  as  good  as  any,  for  the  purpose  that  has 
been  described.  Other  considerations  make  the  conic  spiral  (CS)  trajectory  advantageous 
(see  Chapter  3). 

Since  it  is  evident  that  both  the  PN  and  CLA  trajectories  can  be  taken  as  planar,  the 
following  discussion  is  limited  to  a  demonstration  that  the  planar  CHE  trajectory  is  as  good 
as  any  other  CHE  trajectory: 

Writing  (1.1)  in  the  form 

(r  •  f)^  =  {vr  cos  4>)^  =  (v  cos  (f>)^(r  •  r)  (2.1) 

and  differentiating  with  respect  to  time,  one  obtains 

r  •  r  =  — u^sin^^  (2-2) 

Let  us  find  the  most  general  acceleration  that  determines  a  constant-speed  CHE  trajectory 
as  determined  by  the  conditions 

V  =  Vo,  =  <^o  (2-3) 

Let  the  acceleration  be  expressed  in  the  form 

r  =  Ar  +  Br  +  Cn  (2-4) 

where  A,  B,  and  C  are  constants  to  be  determined,  and  n  is  a  unit  vector  perpendicular  to 
the  plane  of  r  and  r.  It  is  possible  to  express  an  arbitrary  vector  in  the  form  (2.4)  because 
the  base  vectors  used  span  the  three-dimensional  space  except  in  the  trivial  case  where  r  is 
parallel  to  r.  The  conditions  (2.2)  and  r  •  r  =  0  lead  to  the  equations 

r^A  —  (vr  cos  4))B  =  -u^sin^<^  ^2  5^ 

—(vr  cos  <I))A  -f  v^B  —  0 

that  determine  A  and  B,  where  v  and  <j)  are  constants  according  to  the  constant-speed  and 
CHE  hypotheses  (2.3).  The  acceleration  resulting  from  the  values  of  A  and  B  determined 
by  (2.5)  is 

r  =  — r“^[uV  -f-  (vr  cos  ^)i.)]  +  Cn  (2-6) 

This  is  the  most  general  form  for  the  acceleration  that  satisfies  both  the  constant-speed  and 
the  CHE  conditions.  The  ‘constant’  C  has  not  been  determined,  and  it  can  take  any  value 
at  any  time.  The  LA  magnitude  is 
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If  the  optimum  CHE  trajectory  is  defined  as  that  trajectory  which  obeys  both  the  constant- 
speed  and  CHE  conditions  and  requires  the  minimum  LA  at  each  instant,  then  (2.7)  shows 
that  this  trajectory  is  planar  with  the  constant  C  identically  zero. 

In  the  CHE  trajectory,  all  permissible  accelerations  give  the  same  rate  of  approach  to 
the  target,  namely,  r  =  —  u  cos  <j>c.  However,  the  acceleration  that  lies  on  the  plane  of  (r,r) 
has  the  smallest  magnitude  |r|  as  shown  by  (2.7).  In  the  CLA  part  of  the  trajectory,  the 
acceleration  that  lies  on  this  same  plane  gives  the  largest  rate  of  approach  to  the  target. 

Consider  now  the  LA  required  for  the  missile  to  transition  from  its  current  state  to  a 
CLA  trajectory  that  intersects  the  target.  Simple  geometry  shows  that  this  LA  is  given  by 

00X^1=— sin  ^  (2.8) 

r 

which  is  just  twice  the  value  of  |f.|  needed  to  maintain  a  CHE  trajectory,  as  shown  by  (2.7) 
with  C  =  0.  Accordingly,  at  all  times  during  a  CHE  trajectory,  the  LA  being  applied  in  order 
to  maintain  the  CHE  trajectory  is  just  half  that  required  to  transition  to  a  terminal  CLA 
trajectory  that  hits  the  target.  Upon  transition  to  the  CLA  trajectory,  the  LA  is  abruptly 
doubled. 

The  dependence  of  velocity  on  range  in  the  CHE  trajectory  can  be  found  by  expressing 
it  as  a  linear  combination  of  r  and  fc,  the  unit  upward  vector.  The  condition  (1.1)  along  with 
the  requirement  that  the  velocity  magnitude  equal  v  gives 

r  =  ±[r^  —  (fc  •  £)^]'”^^^[— r~^(^  •  sin cos  4>c  (2-9) 

The  trajectory  shape  is  easily  found  in  polar  coordinates  (r,  0)  on  the  flight  plane,  in  terms 
of  which  the  (constant)  speed  is  given  by 

V  =  \/f2  _j-  j~2ff2  (2.10) 

while  the  condition  (1.1)  becomes 

r  =  vcx)S  (f>c  (2.11) 

Equations  (2.10)  and  (2.11)  determine  the  trajectory  in  the  form  6{r),  where 

^(r)  =  00  -  (tan^c)lii(r/ro)  (2.12) 

with  00  representing  the  polar  angle  where  the  range  is  ro.  Some  typical  CHE  trajectories 
are  shown  in  Figures  2-1  and  2-2. 

In  an  actual  engagement,  the  missile  flies  to  the  vicinity  of  the  target  during  the  initial 
flight  phase  during  which  the  HE  generally  decreases  continually.  When  the  HE  becomes 
as  small  as  4>c  there  is  transition  to  a  CHE  trajectory  during  which  the  LA  continually 
increases.  But  before  the  LA  becomes  as  large  as  Oc  there  must  be  another  transition  to  a 
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CLA  trajectory.  This  is  because  there  is  an  abrupt  increase  in  the  LA  upon  this  transition, 
and  it  is  the  increased  LA  during  the  terminal  CLA  trajectory  that  intersects  the  target  that 
must  be  constrained  to  not  exceed  Oc- 

Figure  2-3  illustrates  a  hypothetical  sequence  of  these  events  on  a  plot  of  lateral  accel¬ 
eration  against  HE.  During  PN  both  the  LA  and  the  HE  generally  decrease.  When  the  HE 
reaches  the  value  <^c  there  is  transition  to  CHE,  because  otherwise  the  target  would  not  be 
visible  to  the  sensor.  During  CHE  the  LA  increases  until  half  the  limiting  value  is  reached, 
whereupon  there  is  transition  to  the  terminal  CLA  trajectory.  At  this  transition  the  LA 
abruptly  increases  by  a  factor  of  two. 
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Figure  2-2.  Hypothetical  Evolution  of  Trajectory. 
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CHAPTER  3 

APPROACH  FROM  ABOVE 

Consider  now  the  case  where  it  is  deemed  desirable  to  approach  the  stationary  target 
from  above.  Such  a  quasi-vertical  approach  is  advantageous  when  a  view  of  the  target  is  to 
be  compared  with  an  image  as  viewed  from  above.  A  strictly  vertical  approach  is  impossible. 
Let  an  angnlar  deviation  from  the  vertical,  a,  be  defined.  The  flight  of  the  missile  is  to  be 
constrained  to  lie  within  a  circular  cone,  a  <  ac?  for  specified  Oc-  The  axis  of  this  cone  is 
vertical  and  its  apex  is  located  at  the  target.  A  heuristic  argument  indicates  that  the  best 
trajectory  lies  on  the  cone.  Discussion  of  the  initial  flight  segment  that  brings  the  missile  to 
the  beginning  of  the  “conic  spiral”  (CS)  is  deferred. 

Let  the  motion  be  described  in  terms  of  new  position  coordinates  (/?,  0,  a).  The  relation 
between  these  and  the  Cartesian  coordinates  is 

X  =  p  cos  9 

y  =  psm9  (3.1) 

z  =  p  ctn  a 

The  inverse  relation  is  given  by 

/>  =  Vx2  -1-  j/2 

6  =  arctan(j//z)  (3-2) 

a  =  arctan(Y»2^t-^/2:) 

From  now  on,  o:  is  taken  to  be  constant  and  equal  to  etc.  The  velocity  components  are  given 

by  .  .  •  . 

X  =  pcos  9  —  p9 sin0 

y  =  psm 9  p9  cos  9  (3-3) 

z  =  p  ctn  a  —  pa  csc^a 

The  HE  is  denoted  by  <j)^  which  is  given  by  the  relation 

cos<f>  =  -^  (3.4) 

rv 

Let  us  now  examine  the  motion  governed  by  the  conditions  <f>  =  and  a  =  ac,  where  <^c 
and  ac  are  specified  constants.  Then  (3.4)  gives 

p  =  —V  cos  <j)c  sin  arc  (3.5) 

Using  (3.3)  and  (3.5),  the  velocity  components  can  be  written 

X  =  —V  cos  (j)c  sin  ttc  cos  9  —  z9  tan  ac  sin  9 

y  =  —V  cos  (f>c  sin  ac  sin  9  +  z9  tan  ac  cos  9  (3.6) 

z  =  —V  cos  ^c  cos  ac 
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The  requirement  that  (3.6)  represent  a  constant  speed  v  gives 


from  which  6  becomes 


^2  vhm^<t>c 
z^tan^ac 

(3.7) 

(3.8) 

z  tan  ac 

The  ±  sign  indicates  that  the  missile  can  spiral  in  either  the  clockwise  or  counterclockwise 
direction.  Using  this  result,  the  velocity  components  can  be  expressed  as  a  function  of 
position,  becoming 

X  =  vz~^{—x  cos  cos  ac  T  y  sin  <j>c  ctn  ac) 

y  =  vz~^{—y  cos  <j>c  cos  ojc  ±  X  sin  (f>c  ctn  ac)  (3-9) 

z  =  —V  cos  <f)c  cos  ac 

The  upper  (or  lower)  signs  should  be  used  consistently  throughout.  The  acceleration  com¬ 
ponents  can  be  obtained  as  functions  of  position  by  differentiating  (3.9),  becoming 

X  =  v^z~^{±y  cos  <f>c  sin  <f)c  cos  ac  ctn  ac  —  x  sin^^c  ctn^ ac) 

y  z=  v^z~‘^{=fx  cos  <j)c  sin  <t>c  cos  ac  ctn  ac  —  y  sin^<^c  ctn^  ac)  (3.10) 

z  =  0 


The  acceleration  vector  r  lies  on  the  horizontal  plane  and  is  perpendicular  to  the  velocity 
vector.  Its  magnitude  is 

|r|  =  t;^z“^sin  <f>c  ctn  ac  ycos^^c  sin^ac  -|-  sin^^c  (3-11) 

The  integral  of  the  motion  gives  the  time  history  of  the  trajectory, 

z  =  ZQ  —  vt  cos  (j)c  cos  ac 

e  =  i55i£ln  ( 1  -  (3-12) 

sinac  V  ZQ  J 

This  is  the  CS.  It  appears  that  this  simplistic  and  purely  kinematic  discussion  may  work 
well.  In  particular,  the  constant  v  assumption  may  not  be  overly  far-fetched.  It  is  noted 
that  constant  power  is  being  received  by  the  missile  from  its  uniformly- diminishing  potential 
energy. 


THE  BLIND  FLIGHT  PLANE  (BFP) 

As  before,  it  is  impossible  to  hit  the  target  with  finite  LA  capability.  Therefore,  a 
maximum  LA  value  of  etc  is  specified.  At  some  point,  a  transition  is  made  to  a  constant  LA 
trajectory  (CLA)  that  will  intersect  the  target.  During  this  CLA  trajectory  the  HE  may 
become  less  than  <f)ct  meaning  that  the  missile  flies  blind. 

Let  be  the  polar  angle  at  the  point  of  transition.  The  target  and  the  velocity  at  the 
point  of  transition  lie  on  the  cone,  and  also  on  a  plane  to  be  called  the  BFP,  which  is  tangent 
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to  the  cone  along  the  line  9  —  9\.  Let  be  Cartesian  coordinates  on  the  BFP,  with  the 

r)  direction  being  generally  upward.  The  relations  between  the  BFP  coordinates  and  the 
original  three-dimensional  Cartesian  coordinates  are  given  by 

X  =  T)  sin  a  cos  ^  sin  9i 

y  =  •q  sin  a  sin  ^  cos  9i  (3.13) 

z  =  q  cos  a 

Matching  the  three  Cartesian  components  of  velocity  at  the  transition  point,  indicated  by 
the  subscript  “1”,  gives  the  relations 

jyi  sin  Qc  cos  ^  sin  9i  =  picos  9i  —  pi^isin  9i 

77isin  Qfc  sin^i  +  ^icos  9i  =  pisin^i  +  pi^icos  9i  (3-14) 

7/1  =  pi  CSC  Q!c 

The  acceleration  at  transition  is  not  matched  because  it  can  abruptly  change.  The  terminal 
CLA  trajectory  begins  at  {(,  q)  =  (0,  q\)  and  ends  at  (^,  q)  —  (0, 0).  Matching  the  trajectory 
slope  m  on  the  (^,  q)  plane  at  transition  gives 


qi  picscQc 


(3.15) 


The  geometry  of  the  circular  arc  on  the  (^,  rj)  plane  shows  that  the  initial  rj  must  be  given 

by  2 

qi  =  ^^sin  (f)c  (3.16) 

in  order  that  the  missile  trajectory  hit  the  target,  where  oj  is  the  constant  LA  on  the  BFP. 

On  the  cone  the  flight  is  governed  by  the  relation  <t>  =  <l>c,  while  on  the  BFP  the  flight 
is  governed  by  the  relation  a  =  ai.  According  to  (3.11)  the  LA  on  the  cone  at  transition  is 
given  by 

V  sin  (f) Q  CSC  occ  j  0  1  2  i  *  2  j(  Co  i'7\ 

= - y  cos^ <pc  sm  etc  +  sm  (pc  (3- 1  <  j 

ni  V 


®Tn 


The  HE,  (j>,  during  flight  on  the  BFP  obeys  the  relation 

2u^ 


-cos(^c  —  ^)sin  <f>  =  q 


from  which  it  is  seen  that 


dq  2u2 

—  = - cos{(pc  -  2<f>) 

u<p  ai 


(3.18) 


(3.19) 


Since  q  monotonically  decreases  in  the  BFP,  while  the  right  side  of  (3.19)  is  always  positive  in 
the  appropriate  dornain  of  the  variables,  <f>  continually  decreases. ,  Thus,  the  condition  <i>'><j>c 
is  immediately  violated  upon  transition  to  the  BFP,  and  the  entire  flight  from  transition  is 
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blind.  In  addition,  the  missile  immediately  leaves  the  cone  upon  transition  and  even  if  the 
target  were  to  remain  visible,  it  would  not  be  sufficiently  ‘below’  the  missile. 

The  constraint  on  LA  requires  that  am  <  Oc  s-nd  aj  <  Oc-  Considering  (3.15)  and  (3.16), 
these  conditions  provide  the  following  rule  by  which  the  distance  of  blind  flight,  r}i,  in  the 
BFP  is  determined.  Deflne  the  ratio 


<Zj7l 


-\lcOS^<f)c  + 


f  sin  4>c 
\sin  0!c 


(3.20) 


The  value  =  1  marks  the  transition  from  one  criterion  to  the  other.  Let  a  =  ai  he  the 
value  for  which  //  =  1  for  given  (f>c‘  This  value  is  given  by  the  relation 

sin  (^c  =  tan  oi  (3.21) 


Then 


if  a  >  ai,  ^  <  1?  — ^  o,h  — 

if  a  <  oi,  — »  M  ^  = 


2v^  .  , 

771  =  - sm  (pc 

Oc 

2v^u  . 

7?i  = - sin^c 

ac 


(3.22) 

(3.23) 


The  miss  distance  is  determined  by  the  conditions  at  transition  between  the  CS  and 
the  BF  trajectory,  because  after  the  transition  no  further  measurements  are  possible.  The 
position  r_i  and  velocity  at  transition  are  the  final  values  in  the  CS  trajectory.  The 
subsequent  flight  is  along  a  circular  arc  with  center  at  where  lies  in  the  plane  of  the 
missile  position  and  velocity.  Therefore,  express  in  the  form 


Lc  = 


(3.24) 


where  A.  and  B  are  constants  to  be  determined.  Also,  let  be  the  final  velocity  at  the 
target.  This  final  velocity  is  determined  by  the  conditions  £1  +  tf  =  where  C  is  a 

positive  constant,  and  |£y^|  =  |£il  =  v.  Under  these  conditions,  the  final  velocity  is  given  by 


(3.25) 


The  center  of  curvature  is  determined  by  the  conditions  rc  =  {li  —  r^l  =  b,  ■  tf  =  0,  and 
(£1  —  £^)  •  £1  =  0.  These  conditions  require  that 


£c  =  A£i  + 


(1  -^)£i  -ilii 
— ;;2 


(3.26) 


where  A  is  given  by 


—  2(£i  •  £1)^ 

2[v'^r^  —  (£1  •  £1)^] 


(3.27) 


3-4 


NSWCDD/TN-97/190 


Accordingly,  the  radius  of  curvature  is 

b  = 


vrl 


2y't;2rf-(ri  •fi)2 

while  the  required  acceleration  just  after  transition  is  given  by 

and  its  magnitude  is  _ 


(3.28) 


1“  I 

bl  =  ^ 


(£1  • 


(3.29) 


(3.30) 


If  the  true  missile  state  at  transition  were  known,  its  subsequent  flight  would  be  given  by 

r(t)  =  r j  4-  b<^~^sinwt  +  ^^^”^(1  ~  cosut)  (3.31) 


(see  Groves  et  al,  1994),  where 


(3.32) 


This  trajectory  passes  through  the  origin.  That  is,  under  these  supposed  conditions  of 
knowledge  of  the  true  missile  state  at  transition,  the  missile  hits  the  target.  If  the  true 
missile  state  at  transition  is  known  only  approximately,  the  acceleration  to  be  applied  to  the 
missile  at  each  instant  is  given  by  the  derivative  of  (3.31), 


r{t)  =  —tiUJ  sina;t  +  rjcos  ut 


(3.33) 


This  acceleration  is  computed  from  the  estimated  target  state  at  transition  and  applied  to 
the  true  state,  in  timestepping,  until  the  missile  altitude  reaches  zero,  whereupon  the  miss 
distance  is  given  by  the  horizontal  missile  position  on  the  plane  2  =  0. 

Figure  3-1  shows  a  typical  CS  trajectory  projected  onto  the  horizontal  plane.  Figure  3-2 
shows  the  projection  of  the  same  trajectory  onto  a  vertical  plane. 
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CHAPTER  4 

GUIDANCE  USING  STATE  ESTIMATION 
BY  THE  KALMAN  FILTER 


Tlie  Kalman  Filter  is  a  general  solution  to  the  minimized  mean  square  estimation  problem 
for  a  linear  dynamic  system.  The  Kalman  Filter  gives  a  minimum  error  variance  recursive 
algorithm  to  estimate  an  instantaneous  state  of  a  linear  model  system  by  using  a  measure¬ 
ment  linearly  related  to  the  state,  but  corrupted  by  Gaussian  white  noise  taken  at  discrete 
time  intervals.  For  a  Gaussian  random  variables,  the  Kahnan  Filter  provides  the  minimum 
mean-square  error  estimate  of  the  state  and  for  variables  of  forms  other  than  Gaussian,  the 
Kahnan  Filter  is  the  best  linear  estimator  of  the  state  in  the  mean  square  error  sense. 

To  provide  the  appropriate  acceleration  commands  (3.10)  for  the  missile  to  execute  the 
CS  trajectory,  an  estimate  of  the  missile  position  at  each  instant  is  required.  The  procedure 
is  to  first  apply  a  time  update  to  estimate  the  missile  state  at  a  future  time  given  its  current 
state.  Then  a  measurement  update  is  apphed.  The  time  update  is  usually  based  on  a 
simple  extrapolation  scheme,  such  as  the  assumption  that  the  missile  velocity  will  remain 
constant  during  the  extrapolation  interval.  This  is  the  constant  velocity  (CV)  update.  In 
the  CS  trajectory,  a  time  update  using  a  Constant-Turning-Rate  assumption  will  make  the 
mecisurement  innovations  smaller  and  lead  to  smaller  estimation  errors.  The  model  of  Blair 
et  al  (1997)  is  used,  based  on  the  Constant- Turning-Rate  trajectory  of  Groves  et  al  (1994). 

The  linear  dynamic  model  commonly  used  for  a  target  in  track  is  represented  by 

Ajfc  =  Fk-iXj-^i  +  (4.1) 

where  Xk  represents  the  states  of  the  target  at  time  k.  For  the  constant  Turn-Rate  model, 
the  elements  of  X^  are  position  {x,y,z),  velocity  (x,y,z)  and  acceleration  {x,y,z),  wj;  ~ 
X{0,  Qk)  is  the  process  noise,  assumed  Gaussian,  and  Fk  and  Gk  define  the  dynamics  of  the 
target.  The  linear  measurement  process  is  given  by 

Zk  =  HkXk  +  Vk  (4.2) 

(see,  for  example,  Gelb,  1974)  where  Zk  is  the  discrete-time  measurement  at  time  k,  Vk  ~ 
A^(0,  Rk)  is  the  measurement  noise,  assumed  Gaussian,  and  Hk  is  a  matrix  relating  the  target 
state  to  the  measurements. 

Implementation  of  the  standard  Kalman  Filter  algorithm  requires  that  both  the  states 
and  the  measurement  models  be  linear.  For  the  linear  system  prescribed  by  Equations  (4.1) 
and  (4.2),  the  standard  Kalman  Filter  algorithm  is  given  by 
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L  Time  Update: 

=  Pk-iPk-i\k-iFk-i  +  Gk-iQk-iGk_i  (4.4) 

II.  Measurement  Update: 

Xk\k  =  Xk\k-i  +  Kk(Zk  -  HkXk\k-i)  (4.5) 

Pk\k  =  {I  -  KkHk)Pk\k-i  (4.6) 

Ki  =  fus-iiff  +  Rt)-'  (4.7) 

where  is  the  state  estimate  at  time  k,  Pjt|jfc  is  the  error  covariance  matrix  at  time  k,  and 
Kk  is  the  Kalman  gain  at  time  k.  The  vector  represents  the  predicted  state  estimate 

at  time  k  based  on  measurements  through  time  k  —  1.  and  P;k|jfc_i  represents  the  covariance 

at  time  k  based  on  measurements  through  time  k  —  1. 

The  discrete-time  dynamics  for  constant  turn-rate  model  is  given  by  (4.1)  with 

/Ak  0  0\  /Bk  0  0\ 

Pit  =  0  Ak  0  ,  Gk=[  Q  Bk  0  (4.8) 

\  0  0  Ajfc/  V  0  0  Pi/ 

where 


0 

0  ^ 

(Bk  0 

^it 

0  , 

Gi  =  0  Bk 

0 

AkJ 

Vo  0 

a? 

~^sin(u;T) 

—  cos(a>r) 

) 

cos(u)T) 

a;~^sin(ujr) 

)  - 

uJsin(ajT) 

cos{(jjT) 

/ 

Bk 

II 

CO 

, 

lr2  T) 

rr\ 

Xk  =  {x  X  X  y  y  y  z  z  z}k 

where  T  is  the  time  interval  between  measurements  k  and  -t- 1,  and  co  =  |i:|/|zi|. 


(4.9) 

(4.10) 


Following  Ohlmeyer  (1985),  we  assume  that  the  discrete  measurements  are  available  in  a 
nonhnear  function  of  the  state  Xk,  and  are  represented  by: 


Zk  =  hkXk  -f  Vk 


(4.11) 


where 


hk{Xk) 


(range 
azimuth  angle 
elevation  angle 


/(xl  +  yl  +  zl) 

arctan(yfc/a:fc) 

i^klJ^  +  Vk] 


arctan 


(4.12) 
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where  and  are  the  measurement  errors  at  time  k.  The  measurement  errors  have 

the  property 

E[vk]  =  0,  E[ojv'[]  =  RkSjk 

where  the  measurement  noise  covariance  matrix  is  defined  by 

/<t2  0  0  \ 

fit  =  0  <7?  0  (4.13) 

Vo  0  ajj 

The  process  noise  has  the  properties 

E[wk]  =  0,  E[wjwl]  =  QkSjk 

where  the  process  noise  covariance  matrix  is  defined  by 

/Ef  0  0  \ 

Ot  =  0  El  0  (4.14) 

V  0  0  El) 

where  is  a  diagonal  matrix  whose  elements  represent  the  variance  in  the  Cartesian  co¬ 
ordinates  of  position,  velocity,  and  acceleration,  and  similarly  for  and  These  are 
assumed  to  be  equal,  that  is, 

EI  =  El  =  El=  0  0  (4.15) 

V  0  0  alJ 

Accordingly,  it  is  assumed  that  the  process  noise  in  position  is  not  correlated  with  that 
in  velocity  or  acceleration,  etc.,  and  that  there  is  also  no  noise  correlation  between  the 
coordinates.  The  subscripts  p,  v,  and  a  in  (4.15)  represent  “position”,  “velocity”,  and 
“acceleration” . 


With  the  existence  of  a  nonlinear  measurements  of  the  state  Xk,  we  must  use  the  Extended 
Kalman  Filter  algorithm  which  handles  the  case  of  systems  possessing  nonlinearities.  The 
basic  idea  of  the  Extended  Kalman  Filter  algorithm  is  linearizing  the  nonlinear  measurement 
equations  about  the  most  recent  state  estimate,  and  then  implementing  the  Standard  Kalman 
Filter  algorithm  which  is  given  by  (4.3)  through  (4.7). 


Linearization  is  obtained  by  using  the  first  order  Taylor  series  of  the  nonlinear  measure¬ 
ment  given  in  (4.11)  about  the  predicted  state  estimate  giving 


Zk  —  HkXk  +  Vk, 


Hkij  — 


dHki 

dXki' 


Hk  = 


dhk{Xk) 

dXk 


(4.16) 
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where 


/ 


Hk  = 


<’^k\k-X 
riji-i 
J'ili-l 
,2 


0  0 
0  0 
0  0 


yk\h-l 

^ibtib-l 

r2 

^k\k-l,h 

>‘k\k-iyk\k-l 


0  0 
0  0 
0  0 


H\k-l 

0 

^k\k-l,h 

r2 


0  0^ 
0  0 
0  Oy 


^A:|fc-1  —  +  y\\k-\  +  ^k\k-l  ’  '''k\k-l,k  — 


(4.17) 


4-4 


NSWCDD/TN-97/190 


CHAPTER  5 

CALCULATION  OF  THE  MISS  DISTANCE 


The  miss  distance  (MD)  is  regarded  at  the  appropriate  measure  of  effectiveness  against 
which  the  various  trajectories  discussed  in  previous  chapters  axe  to  be  compared.  These  MDs 
have  been  estimated  only  for  the  conical  spiral  (CS)  trajectories  described  in  chapter  3.  It 
is  assumed  that  the  missile  state  is  imprecisely  known  through  measurements  of  the  target’s 
relative  distance  and  direction  from  the  missile.  The  missile  state  is  estimated  by  means  of 
Kalman  filtering  as  discussed  in  Chapter  4.  The  CS  trajectory  is  obtained  by  applying  an 
appropriate  acceleration  at  each  instant  based  on  the  missile’s  estimated  state,  according 
to  (3.10).  This  acceleration  is  applied  to  the  missile’s  true  state.  During  the  simulation, 
both  the  missile’s  true  and  estimated  states  are  calculated.  Measurements  are  available  only 
during  the  CS  portion  of  the  trajectory;  after  transition  to  Constant  Turning  Rate  no  more 
measurements  are  possible.  The  measurements  are  generated  by  applying  random  errors 
to  the  true  missile  position.  At  the  beginning  of  the  trajectory,  it  is  assumed  that  perfect 
knowledge  of  the  missile  state  is  available.  While  this  assumption  is  not  realistic  and  will 
tend  to  provide  an  overoptimistic  estimate  of  the  MD,  the  effects  will  tend  to  die  out  during 
the  course  of  the  missile  flight  and  with  the  sequence  of  measurement  updates. 


Monte  Carlo  simulations  are  performed  for  each  set  of  parameter  values.  The  standard 
errors  in  range,  bearing,  and  elevation  are  denoted  by  (Tr,u’j,cre.  Each  estimate  of  MD  is 
based  on  100  Monte  Carlo  experiments. 


The  parameters  determining  the  MD  in  each  Monte  Carlo  experiment  are  as  follows: 

4>c  Minimum  allowed  Heading  Error  (HE) 
ac  Half  cone  angle  in  the  CS  trajectory 
V  Missile  speed  (assumed  constant) 

Cc  Maximum  lateral  acceleration  capability 

dt  Time  interval  between  measurements 

(Tr  Standard  error  in  range  measurement 

aj)  Standard  error  in  bearing  measurement 
(Te  Standard  error  in  elevation  measurement 
cTp  Standard  process  error  in  position 
av  Standard  process  error  in  velocity 

(Jo  Standard  process  error  in  acceleration 
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Table  5-1  illustrates  the  effect  of  measurement  rate  and  angular  error  on  miss  distance 
for  the  indicated  parameter  values. 


Table  5-1 


V 

0.5  Mach 

oc  15  g 

(Jr 

5  m 

<7p  3  m^ 

3  m^s”^ 

(tI  3  m^s-^ 

Otc 

15  deg 

10  deg 

Values  of  Miss  Distance  (m) 

dt  — ^  10  ms  50  ms 

80  ms 

0.3  (deg) 

0.3  1.3 

3.6 

0.5  (deg) 

0.5  1.5 

4.4 

1.0  (deg) 

0.9  1.6 

4.8 

2.0  (deg) 

0.8  1-7 

5.5 

3.0  (deg) 

0.8  1.8 

6.1 

4.0  (deg) 

0.7  1.8 

6.5 

10.0  (deg) 

0.6  1.9 

7.3 

Table  5-2  illustrates 

the  effect  of  missile  speed  and 

minimum  heading  error 

Table  5-2 

dc 

10  g 

dt 

10  ms 

<Tr 

5  m 

<^1 

3  m2 

<^b 

15  deg 

3  m2s“2 

O-e 

15  deg 

3  m^s"'^ 

Oic 

15  deg 

Values  of  Miss  Distance 

(m) 

V 

0.4  mach 

0.5  mach  0.8  mach 

<f>c 

10  (deg) 

0.4 

0.9  0.9 

15  (deg) 

0.8 

1.0  1.0 
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Table  5-3  illustrates  the  effect  of  missile  speed  and  minimum  heading  error. 


Table  5-3 


O/Q 

15  g 

dt 

10  ms 

CTt 

5  m 

3  m2 

15  deg 

3  m2s' 

15  deg 

3  m2s- 

Oic 

15  deg 

<f>c 

10  (deg) 
15  (deg) 


Values  of  Miss  Distance  (m) 

V  0.4  mach  0.5  mach  0.8  mach 

0.3  0.5  0.3 

0.3  0.8  0.8 


Table  5-4  illustrates  the  effect  of  missile  speed  and  limiting  lateral  acceleration. 

Table  5-4 


(Tt 

5  m 

dt 

10  ms 

3  m2 

(^b 

15  deg 

3  m2s“' 

15  deg 

<^1 

3  m2s^' 

OCc 

15  deg 

<t>c 

15  deg 

Values  of  Miss  Distance  (m) 


Oc 

10  g 

15  g 

20  g 

V 

0.4  mach 

0.8 

0.4 

0.4 

0.6  mach 

1.1 

0.9 

1.3 

Table  5-5  illustrates  the  effect  of  missile  speed  and  cone  half-angle. 


Table  5-5 


QrC 

15  g 

dt 

10  ms 

CXr 

5  m 

3  m2 

(^b 

15  deg 

3  m2s‘ 

CTe 

<j>c 

15  deg 

15  deg 

<^1 

3  m2s' 

Values  of  Miss  Distance  (m) 


ac  — > 

V 

10  deg 

15  deg 

20  deg 

0.4  mach 

0.5 

0.4 

0.4 

0.6  mach 

0.8 

0.9 

0.9 
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Table  5-6  illustrates  the  effect  of  limiting  lateral  acceleration  and  cone  half-angle. 


Table  5-6 


<f>c 

15  deg 

dt 

10  ms 

(Tt 

5  m 

<^1 

3  m2 

15  deg 

<^l 

3  m^s-^ 

0-e 

15  deg 

3  m^s"^ 

V 

0.5  mach 

Values 

of  Miss  Distance 

(m) 

Uc  — > 

10  g 

15  g 

20  g 

10  (deg) 

0.4 

0.9 

0.4 

15  (deg) 

1.1 

0.8 

0.6 

20  (deg) 

0.7 

1.0 

1.0 

30  (deg) 

0.8 

1.1 

1.3 

Table  5-7  illustrates  the  effect  of  limiting  lateral  acceleration  and  heading  error. 


Table  5-7 


Oic 

15  deg 

dt 

10  ms 

(Jr 

5  m 

3  m2 

15  deg 

3  m2s~ 

(^e 

V 

15  deg 

0.5  mach 

<^1 

3  m2s“' 

(f>c 

10  (deg) 
15  (deg) 


Values  of  Miss  Distance  (m) 


Qc 

10  g 

15  g 

20  g 

0.9 

0.5 

0.4 

0.9 

0.8 

0.6 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 


Results  shown  in  the  previous  chapter  indicate  that  missile  guidance  in  the  CS  trajectory 
can  achieve  satisfactorily  small  noiss  distances  of  less  than  one  meter,  under  appropriate 
conditions.  Reader  is  cautioned  that  these  results  are  based  only  on  kinematical  studies.  If 
this  method  is  to  be  pursued  it  wiU.  be  necessary  to  consider  missile  aerodynamics,  which 
will  tend  to  increase  the  miss  distances. 

The  process  equation  used  was  based  on  the  constant-tuming-rate  model.  This  resulted 
in  considerably  smaller  miss  distances  than  the  constant- velocity  model,  which  was  initially 
tried.  One  might  expect  somewhat  more  improvement  by  modeling  the  process  on  the  actual 
CS  trajectory  being  flown. 

Experiments  with  measurement  errors  of  different  magnitudes  showed  that  the  miss 
distances  are  not  very  sensitive  to  these  errors.  However,  the  miss  distances  are  quite  sensitive 
to  the  measurement  interval.  In  order  for  the  method  to  be  practical,  the  measurement 
interval  should  be  not  much  greater  than  10  ms  to  achieve  miss  distances  of  less  than  one 
meter. 
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